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Abstract 

We present a method to construct a suitable contour deformation in loop momentum space 
for multi-loop integrals. This contour deformation can be used to perform the integration for 
multi-loop integrals numerically. The integration can be performed directly in loop momen- 
tum space without the introduction of Feynman or Schwinger parameters. The method can 
be applied to finite multi-loop integrals and to divergent multi-loop integrals with suitable 
subtraction terms. The algorithm extends techniques from the one-loop case to the multi- 
loop case. Examples at two and three loops are discussed explicitly. 



1 Introduction 



The automation of NNLO calculations is an ambitious project. Among other methods, a purely 
numerical method can be one approach to achieve this goal. To be more specific we think about 
extending the numerical method used at NLO Q]-(9]| towards the NNLO case. The main ingredi- 
ents within this approach would be a set of subtraction terms to render all contributions individ- 
ually finite and a method for the contour deformation such that loop integrals can be performed 
numerically by Monte Carlo methods. For performance reasons this numerical loop integration 
will then be combined with the integration over the phase space of the final particles in a single 
Monte Carlo integration. 

As a first step in this direction we consider in this paper a method for the contour deformation 
of multi-loop integrals. We generalise the method for the contour deformation from the one-loop 
case to the multi-loop case and put special emphasis on the cases of two and three loops, which 
are the most interesting cases for applications. 

As in the one-loop case there are several variants on how the contour deformation can be 
performed. The main difference is given by the fact that some variants use additional Feyn- 
man parameters while others do not. The use of Feynman parameters has the advantage that the 
construction of the deformation is relatively straightforward. This method has already been dis- 
cussed in the literature for the multi-loop case lfTQ] - [T5l . On the other hand, Feynman parametri- 
sation combines a product of n propagators into a single term, which is raised to the power n. 
The power n magnifies the statistical Monte Carlo integration error we get from regions with 
integrable singularities. This can be acceptable for small values of n, but for multi-parton final 
states it is better to develop a method which avoids this problem from the start. This can be 
achieved by deforming the contour directly in loop momentum space without the introduction 
of Feynman parameters. However, in this case the construction of the deformation is more chal- 
lenging. In this paper we present an algorithm which constructs for a multi-loop integral the 
deformation vector in loop momentum space. The algorithm builds upon an algorithm for the 
one-loop case and considers all possible cycles (i.e. one-loop sub-diagrams) within the multi- 
loop diagram. We show that a naive iterative procedure (e.g. picking a one-loop sub-diagram 
and treating the remainder recursively as a ( / — l)-loop diagram) will not work. The sum over all 
cycles is essential. 

The paper is organised as follows: In the next section we introduce our notation and the 
basic concepts. The algorithm for the direct contour deformation in the loop momentum space 
for multi-loop integrals is presented in section [3] Section |4] provides checks and examples. Our 
conclusions are given in section [51 The algorithm for the multi-loop case builds upon a method 
for the one-loop case. A method for the contour deformation in the one-loop case is reviewed in 
appendix lAl 

2 Notation 

In this section we define the notation. We consider a /-loop Feynman graph with m external lines 
and n internal lines. We label the external momenta by pi, p m and the / independent loop 
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Figure 1: The graph of the crossed double box (left) and the associated chain diagram (right). 



momenta by k\, k\. The momenta flowing through the internal lines are denoted by qi with 
1 <i <n. The momenta qi can be expressed as a linear combination of the external momenta pj 
and the loop momenta kj with coefficients —1, or 1: 

/ m 

= L P<i k i - L <*i]Ph py,a y e {-1,0,1}. (l) 

;=i 7=1 

We denote an internal propagator by 

1 

qf — mf + iO 



The infinitesimal imaginary part iO indicates the direction into which the contour should be de- 
formed in the case where the propagator goes on-shell and is not pinched. The integral which we 
want to consider is given by 



It 



d 4 k - 

r R{k u ...M)Y[Dj. (3) 

7=1 



(my 



We assume that the integral is finite. In other words, we assume that the function R(k\,...,kj) 
contains appropriate subtraction terms such that the integral is finite. The function R(ki, ...,kj) 
is either a polynomial in the loop momenta k\, k/, or - more generally - a rational function 
in the loop momenta with poles which are sufficiently far away from the integration contour. 
The latter possibility occurs already at one-loop through the ultraviolet subtraction terms, which 
introduce a new "ultraviolet" propagator with an arbitrary mass ^uv- Beyond one-loop there will 
be more than one ultraviolet propagator. Taking /jy V large enough on the negative imaginary axis 
avoids that the integration contour comes close to this pole. We consider therefore the poles of 
R(ki,...,ki) to be harmless and seek an integration contour, which avoids whenever possible the 
poles of the propagators Dj. The integration contour is entirely determined by the n propagators 



Dj- 



It will be useful to group the internal propagators Dj into chains [Tl6l . Two propagators 
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Figure 2: The quadruple box (left) and two possible chain diagrams (middle and right). The 
propagators labelled 1 and 1' belong to the same chain, as do the propagators labelled 2 and 2'. 

belong to the same chain, if their momenta differ only by a linear combination of the external 
momenta. In the representation of eq. (OQ) two internal propagators belong to the same chain if 
and only if the corresponding matrices p (J are identical up to an overall sign. Obviously, each 
internal line can only belong to one chain. We denote the number of different chains by c. From 
a given Feynman graph we obtain a new graph which we call the chain diagram by deleting all 
external lines and by choosing one propagator for each chain as a representative. We illustrate 
this with a few examples. In fig. (OQ) we show the two-loop graph of the crossed double box and 
the associated chain diagram which is obtained by deleting the four external lines. There are three 
chains, containing three, two and two propagators, respectively. For each chain we only draw 
one propagator as a representative, resulting in the chain diagram shown in fig. (OQ) on the right. 
Since we choose a propagator as a representative for each chain, the associated chain diagram 
for a given Feynman graph is not unique. This is illustrated in fig. © for the quadruple box. In 
the graph of the quadruple box the two propagators labelled 1 and \' belong to the same chain, as 
do the propagators labelled 2 and 2'. If we choose the propagators 1 and 2 as representatives for 
the two chains, we obtain the chain diagram shown in the middle in fig. ©. If on the other hand 
we choose 1 and 2' as representatives, we obtain the chain diagram shown on the right in fig. ([2]). 
Although there can be more than one chain diagram associated to a given Feynman graph, this 
non-uniqueness does not affect our algorithm. It is sufficient to pick one chain diagram out of all 
possible chain diagrams. The introduction of chains can be viewed as a convenient tool to group 
propagators together. The algorithm for the contour deformation will treat all propagators within 
a given chain in the same way. The most general chain diagrams for two-loops and three-loops 
are shown in fig. ©. The two-loop chain diagram consists of three chains, the three-loop chain 
diagram consists of six chains. 

We define a cycle to be a closed circuit in the diagram. We can denote a cycle by specifying 
the chains which belong to the cycle. In the two-loop diagram of fig. © there are three different 
cycles, given by 

C (12) )C (13) )C (23)_ (4) 
Here we used the notation that denotes the cycle consisting of the chains and 
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Figure 3: The generic chain diagrams at two-loop (left) and three-loop (right). The two-loop 
chain diagram consists of three chains, the three-loop chain diagram consists of six chains. 

Similar, there are seven different cycles for the three-loop chain diagram of fig. ©, given by 

£•(123) £(146) £(245) £(356) £(1256) £(1345) £(2346) ^) 

We used the notation that C^ k ' and C^ kl ' denote cycles consisting of the chains C^'\C^\C^ 
and ,C(i\c( k \cV\ respectively. We emphasis the difference between the number of cycles 
of a chain diagram and the number of independent loop momenta of a chain diagram, the former 
being always greater or equal to the latter one. The loop number corresponds to the number of 
independent cycles and not to the number of cycles. This distinction is relevant for the work 
presented here: The algorithm for the contour deformation in the multi-loop case will be based 
on the set of all cycles. We will also show that the restriction to an independent set of cycles will 
in general not be sufficient. 

3 The algorithm for the contour deformation 

We consider a /-loop integral with n internal propagators corresponding to c chains. We have 

l<c<n. (6) 

Internal momenta belonging to the same chain differ only by a linear combination of the external 
momenta. Without loss of generality we can order the propagators Z), and the momenta qi, such 
that 

1. the first / momenta q\ coincide with the / independent loop momenta kf. 

qi = ki, l<i<L (7) 

2. the first c momenta qi correspond to different chains. 

We can associate a loop momentum to each chain. With the ordering as above the loop momenta 
k\, ki are associated to the first I chains. For the remaining (c — /) chains we set 

k i = E P'A' l<i<c (8) 
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Figure 4: The four cycles C^ 123 ), C^ l345 \ C^ 1256 ) and C^ 146 ) of the three-loop chain diagram 
containing the chain 



The contour deformation for a /-loop integral is defined by the deformation of the / independent 
loop momenta 



k r = kr + ikXr (ki, ...Jki) , 1 < r < /, 



(9) 



where all k% is real. Note that K r depends in general on all k\ % and not just k r After this 
deformation our integral equals 



/ri 

J r=\ 



1 d% 



(2kY 



dk" 



dk) 



R{k h ...,k l )Y[D j . 



(10) 



The Jacobian 



dk» 



dk) 



(ID 



is the determinant of a (4/ x 41) matrix and can be computed numerically. The vectors Ki, 
K/ provide the directions for the deformation, while the parameter X determines the scale of the 
deformation. 



3.1 The direction of the deformation 

Let S be the set of all cycles. We denote by S r the sub-set of the cycles which contain the chain 
. In the example of the three-loop chain diagram we obtain for r = 1 the set of the four cycles 

Si = |d 123 ),d 146 ),d 1256 ),d 1345 )} (12) 

This is illustrated in fig. ©. The set S r will be used in the construction of K r : The deformation 
vector K r is given as a sum of deformation vectors K( a ), where each K:( a ) corresponds to one cycle 
in S r : 

k, = J>(«) (13) 
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k(°0 corresponds to an individual cycle. The deformation vector K^ a ) can be obtained with one- 
loop methodes by treating all momenta not belonging to this cycle as external. The construction 
of the deformation vector in the one-loop case is summarised in appendix [A] Eq. (fT3l) reduces 
therefore the problem of the construction of the deformation direction in the multi-loop case to 
the simpler problem of the construction of the deformation in the one-loop case. Eq. (fT3l) is the 
main result of this paper. 

Let us discuss the most important example: For the two-loop chain diagram of fig. © we 
have 



Kl 

K 2 



K (12) +K (13) ? 
K (12) +K (23) 



(14) 



while for the three-loop chain diagram of fig. © we have 



Kl 

K 3 



K (123) + K (146) + K (1256) + K (1345) ? 
K ( 123) + K (245) + K ( 1256) + K (2346) ? 
K ( 123) +K (356) +K ( 1345) +K (2346)_ 



(15) 



The reader might worry that the definition of the contour deformation depends on the choice of 
the independent loop momenta and/or singles out particular cycles (for example the cycle C*- 12 -* 
appearing twice in the two-loop case, while the cycles and C^ 23 ) only appear once). This is 
not the case. If we set in the analogy with eq. © 



K; 



Ep 



I <i<c, 



(16) 



we have in the two-loop case 



and in the three-loop case 





Kl 


= k( 12 )+k( 13 ), 




K 2 


= k( 12 )+k( 23 \ 




K 3 


= k( 13 )-k( 23 ), 


K] 


= K ( 123) 


+ K (146) +K (1256) +K (1345) 5 


K 2 


= K (123) 


+ K (245) +K (1256) +K (2346) 5 


K 3 


= K (123) 


+ K (356) +K (1345) +K (2346) 


K 4 


= K ( 146) 


_ K (245) +K (1345) _ K (2346) 


K 5 


= K (245) 


_ K (356) +K (1256) _ K (1345) 


K 6 


= K (356) 


_ K (146) +K (2346) _ K (1256)^ 



(17) 



(18) 



The deformation in each chain is given by the signed sum of deformation vectors over all cycles 
which contain the given chain. The signs follow from the orientation of the momentum flow. 
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We would like to comment on the sum over the cycles S r in eq. (fT3l) . As an example we 
discuss the two-loop case as in eq. (fTTT) . Suppose that the momentum flowing through one of the 
propagators of the chain becomes soft. Then the deformation vector corresponding to any 
cycle containing the chain goes to zero. In the example at hand this implies 

k( 13 ) -> 0, k^ 23 ) -> 0, (19) 

and as a consequence 

K 3 -> 0, (20) 

as it should be. Suppose now that K^ 12 ) would not be present in eq. (fTTT) . Then also Ki and K2 
would go to zero. However this is not correct. The deformation for the chains and is 
not pinched and K^ 12 ) provides the correct deformation in this case. This shows that deformation 
corresponding to all cycles have to be taken into account. 

Let us consider a deformation vector ic 01 ) corresponding to an individual cycle a. This de- 
formation vector is constructed in such a way that whenever one propagator of the cycle a goes 
on-shell (qf — mf — > 0) we have 2qi ■ > 0. Now pick out any propagator i of the original 
Feynman graph. This propagator belongs to several cycles. Let us denote this set by S r . The 
deformation vector K r for the propagator is the (signed) sum over the deformation vectors i<( a ) 
corresponding to the cycles contained in S r . (The sign takes care of the relative orientation.) For 
simplicity of the argument let us assume that all signs are positive. For each a we have then 
2qi ■ K^ 01 ) > whenever qf — mf — > 0. It follows that this inequality also holds for the sum of all 
K (°0 : We have 2qi ■ K r > whenever qf - mf -> 0. 



3.2 The scaling parameter 

The deformation vectors K r defined above provide the right direction for the contour deforma- 
tion. The deformation is correct, if the deformation vectors would be infinitesimal. However, in 
practice we would like them to be of finite length and in fact to take them as large as possible. We 
can think of starting from an infinitesimal length and scale the deformation vectors to the largest 
allowed length. The largest allowed length is determined by the nearest pole of the propagators 
in this direction. The scaling parameter takes care of this. The definition of the scaling parameter 
is very similar to the one-loop case. The scaling parameter X is given by 



X = 



with 



mm 



l,Al,...,A„,Au V ,...,A 



1*1 ' 



Yj/4 
Yj/4 
Yj/2 



2* J < Yj 
< Yj < 
Yj < 



2Xi 



(21) 



(22) 



and 



Xi = 



and Y,- = 



k- 



m. 



(23) 
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Here we used the notation that the momentum flowing through the j-th propagator is given by 

qj + iKj, (24) 
where both qj and Kj have real entries. The values X\, A, n ensure that we do not hit accidently 

one of the poles of the n propagators. On the other hand, the values ^u V > "•' ^uv ensure that we 
do not hit accidently one of the (harmless) poles of the function R(k\ ,...,ki). We assume that we 
have |5| ultraviolet propagators, one for each cycle. We define the parameters X(j V by 

J 1 : 4Kj |j > Im&iJv) 

Also here we used the notation that the momentum flowing through the j-th. ultraviolet propagator 
is given by 

qj + iXj, (26) 

where both qj and K, have real entries. 



4 Checks and examples 

In this section we test the algorithm for the contour deformation for multi-loop integrals by 
calculating several scalar Feynman integral and by comparing our results with known analytical 
results wherever they are available. We consider scalar multi-loop integrals, where all external 
particles are off-shell. These integrals are infrared finite. In addition we require that the integrals 
are also ultraviolet finite. In particular this excludes diagrams with one-loop self-energy type 
sub-diagrams. Since the examples which we consider are ultraviolet and infrared finite we do 
not need any subtraction terms. The normalisation of the integrals is in accordance with eq. © 




where Dj denotes the propagators as in eq. ©. In the literature there are simple analytical results 
for the two-, three- and four-point ladder diagrams with off-shell external momenta and massless 
internal lines [17] . These diagrams are shown in fig. Q. We denote the /-loop two-point ladder 
integral by B^\p 2 ), the /-loop three-point ladder integral by C^ l \p\,p\,p 2 ) and the /-loop four- 
point ladder integral by D^ l \p\,p\,p\,p\,s,t). We briefly recall the analytical results for these 
integrals. We have 

DV\plpiplpl S ,t) = L^j j*<»(X,Y). (28) 
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p + k< 2 > 





Figure 6: The non-planar two-loop three-point function. 




C< 2> 



Figure 7: The degenerate chain diagram of an three-loop ladder diagram. 



The function <3> can be expressed in terms of polylogarithms 

1 



<*> (/ W) 



» {-\yj\\n 2, -\y/x) 
UX(x,y)fy (j-l)\(2l-j)\ 



Li, 



1 



The definitions of the variables are 



Pi 



X 



X(x,y) 



2 ' 

pi 

2 2 
PlPj 

St ' 
(PI+P2) 



xp(x,y) 



P\ 

pi 

2 2 
P2P4 



Uj(-yp(x,y)) 



• (29) 



Y 
t 



st 



1 — x — y) 2 — 4xy, 



p(x,y) 



(P2+P3) 2 , 
2 

1 —x—y+X 



(30) 



In addition we consider the non-planar two-loop three-point function shown in fig. ©, again with 
off-shell external momenta and massless internal lines. The analytical result for this integral is 
given by JH 



c$\pIpIp 2 3 ) = (c^XpIpIpI)) 2 . 



(3D 



The chain diagram of the three-loop ladder graph is degenerate and shown in fig. ©. As a 
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integral 


numerical result 


analytical result 


fl^GeV" 2 ] 
d 2 )[GeV- 4 ] 
cg } [GeV- 4 ] 
D^GeV" 6 ] 


(-3.59±0.05)- 10- 8 
(-1.80±0.05)-10- u 
(-2.93±0.04)-10- H 
(—5.88 ±0.07) -10- 14 


—3.571 • 10 8 
-1.832- 10 11 
-2.904- 10- 11 
-5.897 -10- 14 


[GeV" 4 ] 
c( 3 )[GeV- 6 ] 
D( 3 )[GeV- 8 ] 


(-7.9±0.5)-10- 14 i 
(-5.3±0.6)-10~ 17 j 
(-7.1 ±0.7) -lO" 19 / 


-8.027- 10- 14 / 
-5.389- 10- 17 / 
-6.744- 10 19 / 



Table 1: Results for the various two- and three-loop ladder diagrams. 



consequence the deformation vectors simplify to 

Kl = k ( 123 )+k( 14 )+k( 125 ), 

K 2 = K (123) +K (245) +K (125) +K (234) 5 

K3 = K (123) +K (35) +K (234) 5 

K4 = K (14)_ K (245)_ K (234) 5 

k 5 = k ( 245 )-k( 35 )+k( 125 ). (32) 

We present the results of the numerical evaluation together with the analytical results for the 
two and three-loop ladder diagrams defined in eq. (|28l and the non-planar two-loop three-point 
integral in table ©. The numerical values of the external momenta are given for the two-point 
functions by 

p = (90,0,0,0) GeV. (33) 

For the three-point functions the external momenta are given by 

Pl = (39.7424,-14.1093,0.102709,20.4908) GeV, 

p 2 = (50.2576,14.1093,-0.102709,-20.4908) GeV, 

P3 = (-90,0,0,0) GeV. (34) 

For the four-point functions the external momenta are given by 

pi = (19.6586,-7.15252,-0.206016,8.96383) GeV, 

p 2 = (26.874,7.04203,-0.0501295,-12.9055) GeV, 

p 3 = (43.4674,0.110491,0.256146,3.9417) GeV, 

p 4 = (-90,0,0,0) GeV. (35) 

In table (OQ) we observe a good agreement between the known analytical results and our numerical 
evaluations. The results for the three-loop integrals have a larger Monte Carlo integration error 
as compared to the two-loop integrals. This is of course expected. 
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(a) (b) (c) 




(d) (e) (f) 



Figure 8: Definition of the various two-loop six-point topologies. 



topology 


our result 


(a) [GeV 


-10] 


(-8.66±0.08)- 


10 


-19 


(b) [GeV 


-lOj 


(-1.17±0.02)- 


10 


-18 


(c) [GeV 


-10] 


(-7.75±0.13)- 


10 


-19 


(d) [GeV 


-10] 


(-1.91 ±0.02) ■ 


10 


-19 


(e) [GeV 


-10] 


(-4.64±0.08)- 


10 


-19 


(f) [GeV 


-10] 


(-1.03 ±0.02)- 


10 


-18 



Table 2: Results for the real part of the various two-loop six-point integrals. 



In addition we include here results on the two-loop six-point functions, again with massless 
internal lines. In this case no analytical results are available. The various two-loop six-point 
topologies are shown in fig. ([8]). The external momenta are defined by 

pi = (12.0588,-1.00017,-2.55373,2.65288) GeV, 

p 2 = (18.6089,-8.9195,6.43508,9.61832) GeV, 

p 3 = (13.8389,-4.73227,-6.55009,-1.55854) GeV, 

p A = (25.5377,20.892,4.32472,-8.89684) GeV, 

p 5 = (19.9556,-6.24009,-1.65597,-1.81582) GeV, 

p 6 = (-90,0,0,0) GeV, 

(36) 

Our results for the real part can be found in table Q. Table © demonstrates that two-loop 
six-point functions can be computed easily. 
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In addition we have checked that the algorithm works as expected by adding masses to the 
internal propagators. 

5 Conclusions 

In this paper we discussed an algorithm for the contour deformation in loop momentum space for 
multi-loop integrals. The method can be applied to finite multi-loop integrals and to divergent 
multi-loop integrals with suitable subtraction terms. The problem is reduced to the one-loop 
case through a sum over all possible cycles. We demonstrated on non-trivial two- and three-loop 
examples that the method works correctly. The techniques developed in this paper will be useful 
for a numerical approach towards NNLO corrections and beyond. 

A The contour deformation for an individual cycle 

In this appendix we discuss the construction of the deformation vector for an individual cycle 
-\y e assume that the cycle consists of propagators. For a cycle we construct the vector 
k( k ), which gives the deformation direction. For the construction of the deformation vector 
we work in a specific Lorentz frame, which we take to be the centre-of-mass frame defined by the 
external momenta of the multi-loop integral under considerations (or more generally the external 
momenta of the scattering process under consideration). We denote the centre-of-mass energy 
by Vs, This is a characteristic scale of the problem. 

An individual cycle can be viewed as a one-loop diagram, where all momenta not belonging 
to the cycle are treated as external momenta. In the following we will drop the super-script a and 
adopt the notation of the one-loop case. We consider the pole structure 

- 1 

I 2 ' ( 3V ) 

j=\ (k-qj) -raj 

with n propagators. The loop momentum is denoted by k, and the origins of the n cones are given 
by the four-vectors q\, q n . The masses are given by m\, m n . The method presented here 
is a minor modification of the one given in ref. [1J. The modifications concern the definition of 
the vectors in eq. (|4Q|) below. In ref. HI we considered the physical situation of a 2 — > (n — 2) 
process. In this situation two external particles (the two initial particles) have a negative energy 
component. In the general situation we have to consider the case of an arbitrary number of 
external particles with negative energy component. In order to simplify the notation we write in 
the following for the real loop momentum k instead of k. The various helper functions occuring 
in the definition of the deformation vector and the default values of the technical parameters are 
summarised at the end of this appendix. 
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A.l The deformation vector 

The deformation vector K is constructed as a sum of two terms and a scaling parameter 

K = ^ de (K int + K ext ). (38) 

K ext is given by 

<Ak) = g,v(c+k v + + c-k v _) (39) 

with 

k± = k-P±. (40) 

and gfj V = diag(l, —1,-1,-1) the metric tensor. The four- vectors P± are chosen such that all 
points q\, q n are in the forward light-cone of P + and in the backward light-cone of P . We 
construct P + iteratively through the following algorithm: We start with the list (q\, ...,q n ) and go 
through the following steps: 

1 . Remove all four- vectors, which are in the forward light-cone of some other four- vector. 
(After this step all four- vectors are separated by space-like distances.) 

2. Find the pair (qi, qj) with the smallest space-like separation, i.e. with the smallest value of 
-(<7/-4/) 2 - 

3. Replace qi and qj by the combined four- vector 

Z+ (qi + qMi-qj) , Z%(x,y) = \ (V+ ^ -gV)) • (41) 

4. Go back to step 1 until only one entry is left. Set the final entry equal to P + . 
The construction of P proceeds analogously: 

1 . Remove all four- vectors, which are in the backward light-cone of some other four- vector. 
(Again, after this step all four-vectors are separated by space-like distances.) 

2. Find the pair (qi, qj) with the smallest space-like separation, i.e. with the smallest value of 

-(qi-qj) 2 - 

3. Replace qi and qj by the combined four- vector 

Z (q, + qj, qi - qj) , (x,y) (s°V " S°V) ) • (42) 

4. Go back to step 1 until only one entry is left. Set the final entry equal to P-. 
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This algorithm constructs P± independent of the order of the list (q\,...,q n ). The four- vector 
(P — P + ) is by construction always time-like. We define 



MP = \J(P--P + ) 2 . (43) 

Hp defines a characteristic scale associated to the cycle under consideration. The coefficients c± 
are defined as 

c± = flh 8T (k„mj,M 2 3 ), (44) 

i=l 

A typical value for the parameter M\ is given by 

M 3 = 0.035 max (/u P , Vf\ , (45) 

where \/I is the centre-of-mass energy defined by the external momenta of the multi-loop integral 
and hp has been defined in eq. (|43~T) . 
K int is given by 



n n 



C = -L^-£c^. + < oft (46) 
'=1 'J=i 

with 

kj = k — qi and ky = k — vy. (47) 
The four-vectors vy are defined by 

1 / (mi — mi) , . \ 

vy = « (,r * , j' <48) 

The coefficients are given by 

n « 

Ci = g (^centre, Yh^DlI^ and C <7 = 8 (^centre, Yl , Mf ) fj^;,/, (49) 

/=1 /=1 

with 

We = ^ (£+ + £-). (50) 
Typical values for the parameters yi and M\ are Yi = 0.7 and 

M 2 = 0.7 max (ji P , \fs j . (51) 
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The factors d^\ are defined by 



d 



U 



h 8+ (k h m 2 ,M 2 
h 8 _(ki,mf,M 2 
max [h§(ki,mj ,M\) ,hQ(-2ki ■ k,M 2 



I = i, mi = 

(<7/-4/) 2 = 0, q Q i<q° v mi = 
(qi~qi) 2 = 0, q°>q% m, = 
otherwise 



The factors djjj are defined by 



dij,i = h e (zij,M 2 ) max[h 8 (k h mj,M 2 ),h e (-2krkij,M 
Zij = (qi-qjf -(mi + mj) 2 . 
A typical value for the parameters M\ is given by 

Mi = 0.035v<f. 

We further have 

K soft = y^a^a: 



(52) 



(53) 



(54) 



(55) 



where the index a sums over some pre-defined directions. We take the pre-defined directions to 
be the four Cartesian coordinate directions: 



K = £ so ft (1,0,0,0), ki =E soft (0,1,0,0), 
k 2 = £ S0 ft (0, 0,1,0), k 3 = E soft (0, 0,0,1), 



(56) 



with Zisoft being an energy scale much smaller then the centre-of-mass scale. A typical value is 
£soft = 0.03 \/$. The coefficients c a are given by 



g(£centre,Yl,M 2 2 ) [ Y[ d al \[ d al 

J=l 1=1 



with 



= max[h 8 (ki,mf,y 2 M 2 ),h Q (2krK a ,y 2 M 2 
d~ t = max [h 8 (ki,mf,y 2 Mj),h Q (-2ki ■ v?,y 2 M 2 )} . 

A typical value for y 2 is y 2 = 0.008. 

It remains to define the scaling parameter X cycle for the cycle. We set 



(57) 



(58) 



Ko 



Kint "T K ex t. 



(59) 
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The scaling parameter X cycle is given by 

^ cle = min[Ui,...AnAcoii], (60) 

with 



Yj/4 
Xj-Yj/2 



2Xj < Yj 

< Yj < 2Xj (61) 
Yj < 



and 



Xj = and Yj = (62) 



Next we define a value A, co n by 



with 



K o / K o 



Kon = ^, (63) 



n n 

i'=l i,j=l a 
i<j 

A.2 Helper functions 

We list here the various helper functions appearing in the definition of the contour deformation. 
The functions h§ + and are defined by 

(±k° - y/P + m 2 ) 

h 8 ±(k,m 2 ,M 2 ) = —± J . (65) 

(±k°-Vk 2 + m 2 ) +M 2 

The function /15 is given by 

(\k°\- VP + m 2 ) 2 
h 8 (k,m 2 ,M 2 ) = —± -J- . (66) 

M*°|- VP + m 2 ) +M 2 



The function hQ(t,M 2 ) is defined by 



*eW = 77^2° (0, (67) 
t + Mf 

where Q(t) is the Heaviside step function. The function g is given by 

^ Mi) = <68) 

The operation "o" denotes the Euclidean scalar product of two four- vectors. 
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